Relative Entropy of a Freely Cooling Granular Gas 
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Abstract. The time evolution and stationary values of the entropy per particle of a homogeneous freely cooling granular 
gas, relative to the maximum entropy consistent with the instantaneous translational and rotational temperatures, is analyzed 
by means of a Sonine approximation involving fourth-degree cumulants. The results show a rich variety of dependencies of 
the relative entropy on time and on the coefficients of normal and tangential restitution, including a peculiar behavior in the 
quasi-smooth limit. 
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INTRODUCTION AND BASIC EQUATIONS 

The dynamics of a dilute granular gas, modeled as a system of hard spheres colliding inelastically with constant 
coefficients of normal (a) and tangential (j3) restitution, can be described at a mesoscopic level by the (inelastic) 
Boltzmann equation for the one-body velocity distribution function /(r,v, <B,f) [1,2, 3,4,5,6,7, 8,9, 10, 11, 12, 13, 
14]: 

^/(r,v,fi),0+v-V/(r,v,fi),0=/[v,fi)|/(r,-,-,r),/(r ) -, •,?)], (1) 
where the collision operator is 

J[v u m\fJ] = o 2 j dv 2 J d(0 2 J dd®(g-a) (g-ff) [(apr^fiVio'OfiVi, (0^-/(^,(0^2,(02)] ■ (2) 

Here, o is the diameter of a sphere, ©(x) is Heaviside's step function, O is a unit vector directed along the centers of 
the two colliding particles, g = Vi — v 2 is the relative velocity of the center of masses, and the double primes denote 

pre-collisional velocities, i.e., (vj' 2 , (0'( 2 ) — \ (vi,2, fl>i a)- Conservation of linear and angular momenta implies 

V1+V2 = V] +v 2 , I(o"-m^-dx\"=I(Oi-m^dx\ u Iwl 2 ' + m^-d x \'{ = IG>2 + m^-d x v 2 , (3) 

where m and / are the mass and the moment of inertia of a sphere, respectively. The pre- and post-collisional relative 
velocities of the two points at contact are g" = g" — jO x ((o" + C0 2 ) an d g = g — §C x + (O2), respectively. The 
coefficients of restitution a and /3 relate the normal and tangential components of g to those of g": 

d-g = -ad-g", dxg = -/3dxg". (4) 

The coefficient a ranges from a = (collisions perfectly inelastic) to a = 1 (collisions perfectly elastic), while the 
coefficient j3 ranges from j3 = — 1 (spheres perfectly smooth) to j3 = 1 (spheres perfectly rough). Except if a = 1 and 
J3 = ±1, energy is dissipated upon collisions. Using Eqs. (3) and (4), it is possible to obtain the (restituting) collision 
rules [8, 11, 12, 13, 14] 

v'/ = vi-C, v 2 ' = v 2 + C, a'{ = flh - — a X C, fl>' 2 ' = fi>2-— ff xC, (5) 

ok 41 ax 

where k = 41 /mo 2 is the reduced moment of inertia and we have called 

C = — - — (a-g)a + — — — - — g-(ff-g)<T--ffx(G)i + a)2 
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(6) 



The basic velocity moments are the number density n(r,t), the flow velocity u(r,f), the average spin (or mean 
angular velocity) £l(r,f), the granular translational temperature r tt (r,f), and the granular rotational temperature 
r rot (r,f). They are defined as 



n{r,t) = jdvjda>f{r,v,a,t), u(r,f) = (v>, Q(r,f) = <»>. T«(r,t) = y([v - u(r,f)] 2 >, T mt (r,t) = '-(co 2 ), 

(7) 

where we have introduced the short-hand notation 



(y/(r,v,<o,f)) = — L- jdv j da>\if(r,v,a>,t)f{r,v,a>,t). 



(8) 



As usually done in the literature on granular gases, the Boltzmann constant is absorbed in the definition of T tr and 
r rot , so that these quantities have dimensions of energy. Conservation of mass and linear momentum imply that 
Jd\J d(QJ[\,G>\f,f] = J dvf d<QvJ[v,a>\f,f] = 0. On the other hand, neither the total angular velocity nor the 
translational or rotational kinetic energies are in general conserved by collisions: 



~dt 



(9) 



coll 



= ljdyj da>e>J[v,a>\f,f] = A, 
dT * ^^JdvJda>v 2 J{v,a>\f,f} = -^T\ ^ =1- Jdvj 'da>co 2 J{v,a>\fJ} = -r t T mt . (10) 



dt 



coll 



The above equations define the spin production vector A and the energy production rates <i; tr and <i; rot . While and 
<i; rot can be positive or negative, the net cooling rate £ = {T lt ^ tv + T mt ^ mt ) / (T a + T mt ) is positive definite [14]. 

A granular gas is intrinsically out of equilibrium and thus, in contrast to the case of energy conservation, the 
evolution of / does not obey an H theorem [15, 16, 17], even if the system is isolated. On the other hand, it is 
worthwhile introducing the Boltzmann entropy density n(r,/)s(r,/) [18], where 



s(r,t) 



In 



/(r,v,<P,f) 
K 



(11) 



is the entropy per particle, K being an irrelevant constant with the same dimensions as /. Although n(r,t)s(r,t) does 
not qualify as a Lyapunov function in the inelastic case, its introduction is justified by information-theory arguments 
[19] and also to make contact with the Boltzmann entropy density of a conventional gas [17]. Given local and 
instantaneous values of the quantities n(r,f), u(r,f), £2(r,f), r tr (r,f), and r rot (r,f), the velocity distribution function 
that maximizes the entropy density is the two-temperature Maxwellian 



/o(v,»)=n 



ml 



3/2 



exp 



m(v-u) 2 /(ffl-Q) 2 



2r tr 



2T 



(12) 



where r rot = (7/3n)((fi) - £i) z > = T mt - m 2 /3. The associated entropy per particle is 



so = 3 — In 



ml 



k V4^ 2 r tr r r ' 



3/2' 



(13) 



where we have taken into account that / dvf d(D /o(v, fi))ln /o(v, fi>) = / dvf d(0 /(v, fi>)ln /o(v, fi)). From Eqs. (9) 
and (10) we obtain the corresponding entropy production as 



ds 
dt 



coll 



§ rot r ot -2/n-A/3 

T mt -IQ. 2 /3 



(14) 



In the case of perfectly smooth spheres (j3 = —1), A = <^ rot = and § tt > 0, so that dso/dt\ coll < due to cooling 
effects [18]. On the other hand, in the more general case of rough spheres (/3 < 1) the Maxwellian entropy production 
dso/dt | coll does not have, in general, a definite sign. 



Obviously, sq > s and thus the amount of missing information associated with the true distribution / is upper 
bounded by Eq. (13). Therefore, the relevant information-theory quantity is the (local and instantaneous) relative (or 
excess) entropy per particle 

As(r,t)=s{r,t)-s {r,t) = -QnR(r,v,o,t)), R(r,y,m,t) = ffi^- . (15) 

/o(r,v,fi),rj 

In principle, the rigorous computation of As(r,t) is not possible unless the solution to the Boltzmann equation (1) 
is known. In order to circumvent this difficulty, one needs to resort to approximations. The aim of this paper is to 
specialize to a homogeneous freely cooling granular gas and evaluate As(t) within an approximate scheme by assuming 
a truncated polynomial expansion for R(y, 03, t ) and keeping terms at the lowest order. As will be seen below, we will 
mainly focus on the impact of roughness (j3 ^ — 1) on the evolution properties of As(t ). 



EVALUATION OF THE RELATIVE ENTROPY 



In general, the ratio R can be expanded in a complete set of orthonormal polynomials {*Pk} as 

«(r,v,»,f) = l+Y\4 k (r,0 , I'k(r,v,fi>,0, 

k 

where k denotes the appropriate set of indices and the orthonormality relation is 

1 



n(r,t) 



d\ / rfO)/o(r,v,fi),f)^(r,v,a>,f)^ k /(r,v,fi),f) = &, 



(16) 



(17) 



This implies that Ak = Thus far, all the equations are formally exact. 

Now, we first assume that / is sufficiently close to fo as to apply the linear approximation (In/?) w (R — 1). In that 
case, Eq. (15) becomes 

As(r,0«-£|A k (r,0| 2 . (18) 

k 

This approximation is consistent with the negative-definite character of the relative entropy A*. Next, we particularize 
to a homogeneous and isotropic system. In that case, u = CI = and the complete set of orthonormal polynomials is 
[20] 



r(fc 1 +2i+f)r(£ 2 +2i + f) 



[Ta)Y(A£+\)h\k 2 \ 



Here, 



c = 



vv = 



03 



m 



y/2T mt /l 



(19) 
(20) 



ice) 

are the reduced translational and angular velocities, while L„ ' (x) and Pu{x) are Laguerre (or Sonine) and Legendre 
polynomials, respectively. Note that ^k^e is a polynomial in velocity of total degree 2(k\ +k2 + 2£). Truncation of 
the expansion (16) after terms of fourth degree yields 



R(v,(D,t) 



l+A 20 o(f) l I , 20o(v,ffl,f)+Ao2o(f) l Po2o(v,(0,f)+Ai 1 o(f) l I , iio(v,(0,f)+Aooi(f) x I , ooi(v,(0,f) 



,020/4 c 2 , 15 \ fl02 / 4 , 2 , !5 

1 + c - 5c + —r + w - 5vv + — 



an\ c - - w - - 



(c-w) 2 



where the cumulants 020, «02> a n, and fo are defined as 

«20 



30 A200 = Ts^ ~ h ° 02 = 7E Am) = T5 {w4) h 



2 4, 2 2 . 2 4 

an = 2 A no = vv ) - 1, = ^=^001 = ^ 



((c-w)Vj(^ 2 ) 



2 2 

~Y~ 

(21) 

(22) 
(23) 



The truncated Sonine expansion (21) implies that Eq. (18) becomes 



MO ~ -y4>(0 - y«m(0 - \a 2 n {t) - 5 -b 2 {t). 



(24) 



Thus, in order to obtain the time evolution of As(t) we need to deal with the evolution equations for the cumulants. 
Taking velocity moments in both sides of the Eq. (1) (with Vf = 0) it is possible to get [14, 20] 



^— = -(£120-^02), 5 T ln(l+a 20 ) = — 5M20-T- 

I i 15 \ 1 +O20 



a T l n (l+a 02 ) = -i('5 jUo2 _^i_ 

ID \ 1 +«02 



, (25) 



<? T ln(l +fln) = - 

y 



3 / x M22 

t(M20 + M02)- T - 

2 1 + a 1 1 



-0*20 + Ate)- — -y- 

2 l+a n + f/? 



• (26) 



Here, dx = Vo{t)dt, where Vo(t) = no 2 yj2T tY (t)/m is a nominal collision frequency, and 

M*,* 2 - -— fdvf d»(*^J[v,a\f,f\, Vb = ~4r / dy / ^(c-w) 2 /[v,fi)|/,/], (27) 
nVo J J nVo J J 

are (reduced) collisional moments. Note that = | Vo JU20 and <i; rot = |vo/io2- 

While Eqs. (25) and (26) are formally exact, they do not constitute a closed set of equations since the collisional mo- 
ments are functionals of the whole distribution function /. On the other hand, in the spirit of the Sonine approximation 
(21), it is possible to insert Eq. (21) into Eq. (27) and neglect terms nonlinear in the cumulants to get 



Hi = C\l } (9) +C^(0)fl2O + C${e)ao2 + C\ l ({e)a n + c£ ,; (0)i, i = 20,02, 11, b 



(0/ 



(28) 



where the twenty coefficients Cj\o) depend on the temperature ratio = T ml /T tr , the coefficients of restitution a 



and j3, and the reduced moment of inertia K. The six coefficients Cq' and Cy } for i = 20,02,/? have been derived in 



-(0 



Refs. [9, 10], while the ten coefficients Q U) and C^, C^, and C\'{ for i = 20,02, 1 1 have been derived in Ref. [14]. 
Finally, the four coefficients and cj^ for j = 20,02, 11 have been derived in Ref. [20]. Insertion of Eq. (28) 
into Eqs. (25) and (26), plus further linearization with respect to the cummulants of jU4o/(l +a2o)> M04/(l +002). 
^22/(1 + fl n). an d + a n + \b), allows us to solve the closed set of equations for 9{t), a2o{t), 1202(f)' a n( f )> an d 
b(t) starting from a given initial condition. Then, the temporal dependence of As(t) is obtained from Eq. (24). 



-(0 r« 



(0 



RESULTS 

We have numerically solved the set of equations (25) and (26) with the Sonine approximation (28) for several values 
of a and j3. In all the cases the spheres are assumed to have a uniform mass distribution (i.e., K = |) and the initial 
condition corresponds to 0(0) = r rot (0)/r tr (0) = \, a 2 o(0) = 002(0) = a n (0) = b(0) = 0, so that As(0) = 0. 

Figure 1 shows — As(t) for a representative number of pairs (ot,/3). As is commonly done, the temporal evolution 
is monitored not by time t but by the accumulated number of collisions per particle, # collisions/particle = ^/2kt. A 
log-log plot is used in Fig. 1 because of the great disparity in the orders of magnitude of both \ A(s) | and the relaxation 
time for the different pairs (a,j3). For perfectly rough spheres (j3 = 1), As exhibits a non-monotonic behavior and 
reaches its stationary value As(°°) after less than 10 collisions/particle. Note that in the conservative case of perfectly 
elastic and rough spheres, (a,j3) = (1, 1), As(°°) = 0. For j8 = the relaxation is slightly slower and the stationary 
value |As(°°)| is clearly higher than for j3 = 1. As the spheres become less rough, a larger number of collisions is 
needed to activate the rotational degrees of freedom and, therefore, the relaxation period toward the stationary value 
becomes longer. In fact, as shown in Fig. 1, the stationary value corresponding to j3 = —0.99 has not been reached 
after 10 3 collisions/particle. We have observed that, in general, the relative entropy is dominated by the contribution 
~TT fl 02( f ) °f Eq. (24). This implies that the deviation of s(t) from so(t) is mainly due to the kurtosis of the angular 
velocity distribution function, the other cumulants playing a smaller role. 

Figures 1(a) and 1(b) also include the evolution of — As(t) for perfectly smooth spheres (j3 = —1). In that case, 
the angular velocity of each individual particle is not affected by the dynamics and thus it does not play a role 




0.1 1 10 100 1000 ' 0.1 1 10 100 1000 0.1 1 10 100 1000 



# collisions/particle # collisions/particle # collisions/particle 

FIGURE 1. Log-log plot of — As(t) versus the number of collisions per particle for (a) a = 0.8, (b) a = 0.9, and (c) a = 1. The 
curves in each panel correspond to j8 = 1, 0, —0.9, and —0.99. Moreover, the special case of perfectly smooth spheres (j3 = — 1) is 
included in panels (a) and (b). The dashed horizontal lines represent the stationary values. 
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FIGURE 2. Plot of the stationary values -As(oo) versus j8 for a = 0.8, 0.9, and 1. 




different from that of a label or tag. Therefore, if J3 = — 1 the rotational degrees of freedom are irrelevant and Eq. 
(24) must be replaced by As sm (f) w — -^a^oM- As can ^ e seen f rom Figs. 1(a) and 1(b), the evolution of As(t) for 
P = —0.99 is practically indistinguishable from that ofAs sm (f ) for a transient period longer than the relaxation time 
of the latter quantity. However, as time progresses, the temperature ratio r rot /r tt becomes so large that the rotational 
and translational degrees of freedom become strongly coupled and As(t) for j3 = —0.99 eventually departs from 
As sm (°°). A similar effect has been previously reported for binary mixtures [13]. It is important to remark that even the 
translational cumulant fl2o(°°) differs for smooth spheres (/3 = — 1) from the one obtained for quasi-smooth spheres in 
the limit /3 -> -1 [14]. Note that As sm (f) = if (a,j3) = (1,1) and so it is absent in Fig. 1(c). 

The stationary values — As(°°) are plotted versus j3 for a — 0.8, 0.9, and 1 in Fig. 2. As can be seen, |As(°o)| does 
not present a monotonic behavior with respect to j3 and tends to increase with decreasing a. It is also apparent that a 
rapid change of behavior takes place in the quasi-smooth region J3 ~ — 1 . 



CONCLUSION 



In this paper we have studied the time evolution and stationary values of the Boltzmann entropy of a homogeneous 
freely cooling granular gas, relative to the maximum entropy consistent with the actual translational and rotational 
temperatures. Despite the lack of an H theorem for granular gases, the quantity As(t) measures the amount of missing 
information on the microscopic state of the gas and, consequently, represents an insightful tool to characterize the 
departure of the velocity distribution from a two-temperature Maxwellian. The evaluation of As(t) has been carried 
out by expanding the velocity distribution function in orthogonal polynomials, truncating after fourth degree, and 
keeping the lowest order in the coefficients. As shown by Figs. 1 and 2, the results show a rich variety of dependencies 
of As(t) on time, a, and j3, including a singular behavior in the quasi-smooth limit J3 — > — 1. We plan to undertake in 
the near future a comparison of the results derived in this paper with computer simulations. 
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